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Abstract 

Double tunneling junctions of ferromagnet-superconductor-ferromagnet 
electrodes (FSF) show a jump in the conductance when a parallel mag- 
netic field reverses the magnetization of one of the ferromagnetic elec- 
trodes. This change is generally attributed to the spin-valve effect or to 
pair breaking in the superconductor because of spin accumulation. In 
this paper it is shown that the Meservey-Tedrov effect causes a similar 
change in the conductance since the magnetic field changes the energy 
spectrum of the quasi-particles in the superconductor. A reversal of the 
bias reverses the sign in the conductance jump. 

PACS: 73.50.-h, 73.50.Bk, 73.23.-b, 73.25.+i 

1 Introduction 

During the last five years single electron transistors (SET) with two ferromag- 
netic electrodes and a superconducting island have been studied experimen- 
tally P, [2], El, 0] and theoretically UJ|. Experimentally one generally uses 
ferromagnetic-superconductor-ferromagnetic double junctions which consist of 
an Al strip of length of about Ifim, width of 50-100nm and thickness of about 
20nm. The Al is oxidized, and two Co electrodes with slightly different width 
and twice the thickness cross the Al strip at a separation of a few lOOnm. They 
form two tunneling junctions. Fig.l shows the schematic arrangement of the two 
ferromagnetic Co electrodes and Al island. A magnetic field is applied parallel 
to the Co strips and aligns the magnetization of the two Co strips. Then the 
magnetic field is reversed. At a magnetic field B sw the wider Co strip flips its 
magnetization to be parallel to the magnetic field while the narrower Co strip 
remains anti-parallel to the external field because its coercitive field is larger. 



At the same time the current through the double junction changes abruptly at 
B sw . At the field B sn the narrower Co strip also reverses its magnetization and 
the magnetizations of the two Co strips are again parallel to the external field. 
(The relative orientations of the magnetic field and the magnetization of the two 
Co electrodes is shown in Fig. 4). If one applies constant bias to the junction 
then the current shows a jump at each of the fields B sw and B sn (with opposite 
sign). Such jumps at the fields B sw ,B sn have been observed in a number of 
experiments pQ, [3], @- 

In the theoretical discussion one generally considers two mechanisms which 
change the current (i.e. conductance) of the double junction in the field range 

sw 5 ^sn ) ■ 

1. Spin-valve effect: When the magnetizations of the Co strips ni! and m 2 
are both parallel to y then one has a large density of states in both Co 
electrodes for the spin moment up electrons, while the spin moment down 
electrons have a small density of states in both electrodes. For the (spin) 
moment up one has two small resistances Rt] in series and for the other 
direction two large resistances Rt[. The total conductance is then G] \ — 
1/ (2i?tf) + 1/ (2Rt[). If the two Co strips have opposite magnetization 
then the conductance is G-| \ = 2/ (Rf\ + Rti)- It is easy to show that 
G] \ > G-\ [. Therefore the current should drop inside the field window 

2. Gap reduction due to spin accumulation: In the anti- ferromagnetic align- 
ment one obtains spin moment accumulation because the spin moment 
up electrons have a small resistance for tunneling onto the island and a 
large resistance to tunnel off the island while the opposite is true for spin 
moment down electrons. This spin moment accumulation can reduce the 
superconducting gap of the Al island. This will lead to an increase of the 
conductance in the field window (B SW1 B sn ). 
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Fig.l: The schematic geometry of an FSF double junction, 
consisting of Co/Al/Co. 

In this paper we want to show theoretically that there is an additional con- 
tribution to the current because of the Zccman effect which shifts the excitation 
spectrum of the quasi particles in the Al by /i e B (/i e =moment of the spin up 
and down electrons, B=external magnetic field). This effect has been inten- 
sively studied by Meservey and Tedrow in many beautiful experiments (see the 
review article [§])• In a series of papers ^Oh [TTj - [T^j their group investigated 
the tunneling I-V-curve for ferromagnet-superconductor tunneling junctions in 
different magnetic field. They showed that the I-V-curves were asymmetric with 
respect to the voltage (because of the different density of spin up and down elec- 
tron at the Fermi surface). From the asymmetry they derived the polarization 
of the effective density of states of the tunneling electrons. To our knowledge 
the magnetic field and the magnetization were always parallel to each other in 
their measurements. 

In this paper we want to demonstrate that one has to take the field and 
spin dependence of the quasi-particles in the Al into account when calculating 
the current through the double tunneling junction. To demonstrate this effect 
we consider a double tunneling junction in which the two junctions are so far 
separated that the spin-orbit scattering destroys any spin polarization along the 
diffusion path of the electrons from the first tunneling junction to the second 
one. This means that only the total current 1\ through junction one must 
be equal to the total current I2 through junction two; the spin up and down 
currents through the two junctions can be quite different. 
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2 Theory and Simulation 



2.1 Single junctions 

We first consider a single ferromagnet-superconductor tunneling junction at zero 
temperature. In Fig. 2 the density of states for both metals is plotted after lifting 
the energy bands of the ferromagnet by eU. 



JlLt 

^=^B 



B 
1 



1 

B 



Mt7 



B 
1 



Fig. 2: The tunneling density of spin moment up electrons in 
an FS tunneling junction for different orientations of the 
magnetic field and the magnetization m of the ferromagnet. 

a) B and m are parallel, both pointing in y, 

b) B and m are anti-parallel, B pointing in — y, m in y 
c,d) The bias is reversed. 

In large body of experiments Merservey and Tedrow [HJ showed that a mag- 
netic field parallel to the tunneling junction shifts the excitation spectra of spin 
up and down electrons in the superconductor by /i e B in opposite directions. 
This enhances the current of the majority spin (see Fig. 2a) when the electrons 
are flowing from the ferromagnet to the superconductor. One obtains a spin 
current (with moment up). The I-V-curve is not (point-) symmetric about the 
origin. 

The tunneling current is for spin moment up and down is given at zero 
temperature by 



J T = CN M N ; 



CN m N s 



eU 
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Here Nm and N m are the majority and minority density of states in the 
ferromagnet and N$ is the density of states of the superconductor in the normal 
state. In the superconducting state in the presence of a magnetic field Ns is 



modified by the factor (E + jf e B) /J (E + /~^B) 2 - A (B) 2 . The constant C 



contains the tunneling matrix elements and universal constants. The energy gap 
is given by A (B) . For thin films and stripes which are aligned parallel to an 
external magnetic field we use the result from reference [T^ for the dependence 
of A on the magnetic field : 



where the field B c (T) is determined by the ratio of the penetration depth A (T) 
and the film thickness. 



B c i, (T) is the thermodynamic critical field. 

The use of the density of states in the tunneling current is a dramatic over- 
simplification since the tunneling probability of electrons at different parts of the 
Fermi surface depends strongly on the direction of their group velocity relative 
to the tunneling barrier. So Nm , N m and Ns have to be interpreted as " effective 
tunneling densities of states" . In the present paper we only need the relative 
magnitudes of Nm and N m which are given by the experimental polarization of 
the tunneling current. 

Merservey and Tedrow obtained a number of interesting results for an FS 
junction in a parallel magnetic field: 

• The I-V-tunneling curve is not (point) symmetric about the origin. 

• The tunneling current is polarized and the polarization can be evaluated. 

• The polarization is always parallel to the majority moment of the fer- 
romagnet and not proportional to the d-density of states at the Fermi 
surface. 

They obtained a polarization of 35% for Co/Al junctions. 

There is another interesting consequence of the energy shift of the Zeeman 
effect. Let us consider a single Co/Al tunneling junction, i.e. the left half of 
the Fig.l. We align the magnetic field parallel to the Co strip in the negative 
y-direction and keep the voltage across the junction constant. For simplicity 
we assume that the temperature is (close to) zero. We start with the magnetic 
field — B c (0) which suppresses superconductivity in the Al completely. Then 
we sweep the magnetic field towards +B C (0). As soon as the magnetic field 
takes positive values the magnetization of the Co and the field are anti-parallel 
and therefore the junction is in an instable energetic state. Due to its coercitive 





B c (T) = V^A^r 1 Bo (T) 
a 
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field B sw the Co can maintain the anti-parallel orientation up to the field B aw . 
Then the Co film will switch its magnetization. As a consequence the tunneling 
current will also change. 

In Fig. 3 we calculate the current through the junction using the following 
parameters: the energy gap in Al at zero temperature A = 0.2meV, the field 
that suppresses superconductivity completely B c = 1.5T, the switching field 
B S w = 0.16T, the polarization of the effective density in the Co p — .35. We 
sweep the external magnetic field from — 0.5T to 0.5T. When the magnetic 
field changes sign the Zeeman term changes sign as well. At the magnetic field 
B sw the direction of the Co magnetization m becomes aligned parallel to the 
magnetic field. At the same time the current jumps to a higher value. 

eU (meV) 




B(T) 

Fig. 3: The simulated current through an FS-tunneling junction 
while the magnetic field sweeps from — 0.5T to +0.5T. 
The numbers at the curves give the different bias. At 
B = 0.16T the magnetization of the Co strip reverses. 

In Fig. 3 the calculated tunneling current through a Co/Al junction is plotted 
for constant bias as a function of the magnetic field. The different curves are 
for different bias which is given in meV at the right side of the curves. One 
recognizes that the current shows a jump at the switching field B sw — 0.16T. 
For positive bias the current increases at the switching field while for negative 
bias the (absolute value of the) current decreases. Furthermore the minimum 
of the the I-B-curve is not at B = but shifted to positive field values. If the 
field is then swept from B m = 0.5T to — 0.5T the resulting current curves are 
just a mirror image of the shown curves. 

It is important to note that a reversal of the applied voltage corresponds 
to a tunneling of electrons from the superconductor to the fcrro-magnet (see 
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Fig. 2b, c). In this case the current is smaller if m is parallel to B because for 
an electron with moment up to tunnel from S to F a Cooper pair has to split 
and the moment down electron is elevated by A + \lbB into an excited state in 
the superconductor while the moment up electron tunnels into the ferromagnet. 
The contribution of moment up electrons to the tunneling current is reduced to 

CNfrfNgJ {{eU — [IbB)) 2 — A (B) 2 . Therefore the current changes to a smaller 
value when the magnetization flips. 

2.2 Double junctions 

We now consider the double junction in Fig.l. However, we use a long Al strip 
so that the two junctions are relatively far apart. As we discussed above the 
tunneling current through a single junction Co/Al is polarized. This means 
that one has a source of polarized electrons in the Al strip. Meakawa et al. 
8 calculated an (opposite) shift in the chemical potential for (spin) moments 
up and down. The lifetime of an electron in a given spin state is limited by 
the spin-orbit scattering and the spin state decays as exp (— t/r s f) where T s f is 
the spin- flip lifetime which is 3/2 times the spin-orbit scattering time t so . The 
injected spin polarization at the left tunneling junction diffuses through the Al 
strip towards the right junction and vice versa. Along this diffusion path the 
spin decays. For the resulting spin diffusion length one finds different orders of 
magnitude in the literature, varying between 10-100nm and |14| . 

As shown in Fig.l the total potential difference between the right and the 
left electrode is 2eU. We consider the bias as positive when the potential on the 
right electrode is positive. Then the electrons flow from the left to the right side 
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as shown in Fig. 4. 
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Fig. 4: The current of spin moment up electrons through an FSF-double junction. 

a) The moments mi, m.2 and B are all parallel, pointing in +y direction, 

b) the magnetic field has changed to the — y direction, 

c) the moment m 2 has switched at B sw to the — y direction, 

d) the moment mi has switched at B sn to the — y direction, 
mi, \xi2 and B are all parallel, pointing in — y direction 

In contrast to the goal of the spin valve we choose a large separation of the 
two junctions (in our virtual experiment) so that each junction is unaffected by 
the polarized current of the other. In this case we have excluded all spin-valve 
effects in our virtual experiment. Mathematically this requirement is expressed 
by the condition that only the total currents through junctions Ji and J 2 must 
be identical; the individual spin currents can be different. As a result both spin 
directions experience the same shift in the chemical potential. 

We consider first the special case (a) in Fig. 4 where mi, m 2 and B are 
all pointing in the positive y-direction. In general the currents through the 
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junctions Ji and A are not identical when their bias is the same, i.e. eU. 
Therefore the chemical potential of the island will shift by <j> (which has to 
be determined self consistently). Then the (spin) moment up current through 
junctions Ji and A are given by 



/ 1T =CN s N m ^j({V e + (jy + ^BB) 2 - A 2 ) (1) 
At - CN s N m J ((K - (f> - ii B B) 2 - A 2 ) 

(The symbol f stands again for spin moment up.) 

The other current contributions can be obtained from these currents by 
applying simple rules: 

1 . The contribution of spin moment down electrons is obtained changing the 
sign of the term \ibB in A and A and exchanging Nm and N m . 

2. If mi points in the —y-direction one has to replace Nm by N m in A- 

3. If 1112 points in the —y-direction one has to replace Nm by N m in A- 

4. If B points in the — y-dircction one has to change the sign of the term 
(IbB in A and A- 

It is sufficient to calculate the current for spin moment up (equation £3)) in 
the alignment of Fig. 4a. Then the above rules yield the current for moment up 
and down under all circumstances. For example the corresponding spin moment 
down currents are 



Ai = CN s N M] J((V e + 0-ti B B) 2 (2) 
hi = CNsNmJ ((V e -<P + fi B B) 2 - A 2 ) 

We calculate the total current perturbatively. For sufficiently large bias, 
i.e., eU > A the terms ecj) and [IbB are small compared to y/V^ — A 2 and we 
can expand the different current contributions as a Taylor series in terms of ecj) 
and \ibB up to second order. Since the current depends on the orientation of 
three vectors, mi, m 2 and B we choose the y-direction as a reference direction. 
The value of B is negative when B is anti-parallel to y. We calculate (in 
perturbation) the currents I] [, Jj \, Jj [ for the four possible orientations of 
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mi and m 2 . The results are collected in the following equations 



(3) 



Here the indices of the currents give the direction of the magnetizations mi 
and m 2 with respect to the y direction, for example J| f is the current for mi 
anti-parallel and m 2 parallel to y. When mi and m 2 are anti-parallel to each 
other (i.e. for 1^ ^ and J| f) then the chemical potential of the island e<j) is zero. 
In the parallel orientation one obtains 



e<j> = fOT J 1 t 

e^+^{feS for J n 



The dependence of the currents on the quadratic term ([i B B) 2 is rather 
weak, so that it is sufficient for a qualitative discussion to restrict ourselves to 
the linear dependence on (fisB). When we start with a large negative magnetic 
field (B 1 [ y) then both magnetizations, mi and m 2 , are anti-parallel to y. In 
the linear approximation the current is I\ \ « \/Ve A 2 (Nm + N m ). At the 
positive field B sw the magnetization of the junction J 2 flips and aligns parallel to 
the field. Then the new current is 7j \ which corresponds to a relative decrease 
of the current 

AJ (N M ~ N m ) Ve 
I (N M + N m ) (K 2 - &f B 

At the higher field B sn the other electrode Ji also aligns parallel to y. Then 
the current returns to the original curve since 1^ \ = I\ |. in this approximation. 
It is important to point out that the jump in the current is positive when the 
magnetization of junction Ji (i.e. mi) flips first. Then the current changes from 

In Fig. 5 we sweep the magnetic field from — 0.5T to +0.5T. The current 
through the double junction is plotted versus the sweeping magnetic field. The 
junction J 2 has the switching field B sw = 0.14T while junction Ji has the larger 
switching field of B sn = 0.18T. We use different bias voltages in the range of 
—0.2meV < —ell < 0.7meV. The I — i?-curve in Fig. 5 shows a downward 
displacement in the field window (B sw , B an ). For negative — elf (i.e. U > 0) 
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the absolute value of the current increases. This means that the displacement 
changes sign when the bias is reversed. 




Fig. 5: The current through an FSF-doublc junction while the magnetic field 
sweeps from — 0.5T to 0.5T. The numbers at the curves give the different 
bias. The switching fields for the two Co strips are 0.14T and 0.18T. 

3 The single electron transistor 

When the size of the two tunneling junctions is in the nano-meter scale then 
their capacity is small and the tunneling electrons change the Coulomb energy 
on the island US], C3, [IBI, CU, ED] If an electron from the left electrode 
with the band energy el tunnels onto a state on the island with the band energy 
El and changes the number of electrons on the island from n to (n + 1) then 
conservation of energy requires that 

e 2 Cr> 

E L + Ue = Ei + (2n + 1) — - -T^eU G 

On the other hand if an electron from the island with the band energy £/ 
tunnels into a state on the right electrode with the band energy er and changes 
the number of electrons from n to (n — 1) one has to fulfill the condition 

e 2 C G 

ei + Ue = e R - (2n - 1) — -preUG 

Here Uq is the gate voltage and C^ 1 = Cf 1 +C2 1 +Cq X where C\, C2 and Cq 
are the capacities of the two tunneling junctions and the gate. The Coulomb 
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blockade energy is given by Ecb — e 2 /2Cs- In the following we consider only 
zero gate voltage. 

At a sufficiently large bias (2eU > 2 (A + Ecb)) the island can gain or loose 
up to n electrons where (at T — 0) n is given by n =int [i ((ell — A) / E cb — 1)] 
The probability for n excess electrons on the island may be p (n) which will be 
determined self-consistently. 

First we calculate the currents for spin moment up and mi and ni2 parallel 
to y. For tunneling from the left Co electrode onto the Al island with n excess 
electrons (prior to the tunneling) the current (n) of moment up electrons is 

7i T (n) = p (n) CN M N s \j((eU - (2n + 1) E cb + hbB)) 2 - A 2 

The current from the Al island with n excess electrons onto the right Co elec- 
trode is 

7 2T (n) = p (n) CNmNs^J ((eU + (2n - 1) E Ch - n B B)f - A 2 

The occupation probabilities p (n) are obtained by the condition that the flow 
of electrons on the island with n excess electrons is equal to the out-flow (see for 
example the review article |18j'). This yields simple linear equations for p(n). 
The currents for spin moment down and different orientations of B, mi and m 2 
are obtained by applying the rules which we stated above. The results of this 
calculation are plotted in Fig. 6. We use for the Coulomb energy Ecb the value 
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E Cb = O.lOlmeV. 
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Fig. 6: The current through the FSF-single electron transistor 
for different bias. The switching fields are identical to Fig. 5, 
field sweep is from +1.5T to — 1.5T. Energy gap and Coulomb 
energy are 0.2meV and O.lmeV. 

There are a few kinks in the current curves of Fig. 6 as a function of the 
magnetic field. They occur when when the maximum number of electrons on 
the island changes by one. The magnetic field lowers one subband of the super- 
conductor and reduces the energy gap. Whenever \{eU + /Ug-B) — A (B)] /Ecb 
crosses an odd integer (2n + 1) as a function of increasing \B\ the maximum 
number of electrons on the island increases by one. Furthermore one observes 
that again for negative bias the sign of the relative current jump in the window 
(B sw ,B sn ) changes sign. 



4 Discussion and Conclusion 

In the discussion of a single ferromagnet-superconductor junction we arrived at 
the following conclusions. 

• For electron flow from the ferromagnet F into the superconductor S the 
current increases when the magnetization m aligns parallel to the mag- 
netic field. 

• For electron flow from the superconductor S into the ferromagnet F the 
current decreases when the magnetization m aligns parallel to the mag- 
netic field. 
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From these facts it follows that the current jump in a double junction changes 
sign when one reverses the bias. When the source electrode (the electrode 
from which the electrons tunnel into the island) flips its magnetization first 
then the conductance of the source-island junction increases and therefore the 
current through the SET increases. When the drain electrode (the electrode into 
which the electrons tunnel from the island) flips its magnetization first then the 
conductance of the island-drain junction decreases and therefore the current 
through the SET decreases. Since a reversal of the bias exchanges source and 
drain one finds that the relative change of the current at the field B sw has the 
opposite sign. 

In a nut shell: a flip of the magnetization in the electron source yields an 
increase of the current and a flip of the magnetization in the electron drain a 
decrease. The / — -B-curves for opposite directions of the magnetic field sweep 
are mirror images of each other. 

In this paper we intentionally excluded a spin-coupling between the two 
tunneling junctions. Such a coupling has been observed for example in the 
beautifull spin precession experiment by Jedema et al. 0]. The Meservey- 
Tedrow effect is an additional phenomenon which has to be included in the 
analysis of FSF-single electron transistors. It can be distinguished from the 
effect of spin- accumulation and gap reduction because it changes sign when the 
bias voltage is reversed. 

Abbreviations: SET=single electron transistor, FSF=ferromaget-superconductor- 
ferromagnet. 
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